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Indirect exchange of magnetic impurities in zigzag graphene ribbon 
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We use quantum Monte Carlo method to study the indirect coupling between two magnetic impurities on the 
zigzag edge of graphene ribbon, with respect to the chemical potential fi. We find that the spin-spin correlation 
between two adatoms located on the nearest sites in the zigzag edge are drastically suppressed around the zero- 
energy. As we switch the system away from half-filling, the antifcrromagnctic correlation is first enhanced and 
then decreased. If the two adatoms are adsorbed on the sites belonging to the same sublattice, we find similar 
behavior of spin-spin correlation except for a crossover from ferromagnetic to antifcrromagentic correlation in 
the vicinity of zero-energy. We also calculated the weight of different components of d-electron wave function 
and local magnet moment for various values of parameters, and all the results are consistent with those of 
spin-spin correlation between two magnetic impurities. 

PACS numbers: 73.22.Pr, 75.30.Hx 
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I. INTRODUCTION 

Zigzag graphene ribbon has attracted enormous atten- 
tion due to its potential application in nanoscale spin- 
tronics devices. 1-3 It has zero-energy localized states at 
Fermi level, while those with armchair edges have no such 
peculiar states 4 ' 5 . 

Due to the edge states in zigzag graphene ribbon, the 
magnetic property of the Anderson impurities would be 
greatly modified and shows novel behavior differing from 
that in bulk graphene and in normal metals. When 
two magnetic impurities are located sufficiently far away 
from each other, they can be modeled by a single im- 
purity Hamiltonian 6 ' . However, if they are close to- 
gether, the "spin compensation clouds" of surrounding 
conduction electrons will overlap, and interesting effects 
show up due to the indirect Ruderman-Kitter-Kasuya- 
Yosida (RKKY) coupling 8-10 between magnetic impuri- 
ties, which is mediated by the conduction electrons. 

In the bipartite lattice, the oscillatory RKKY interac- 
tion becomes commensurate and the effective exchange 
coupling is always ferromagnetic (FM) for the same sub- 
lattice but antifcrromagnctic (AFM) for opposite sub- 
lattice at the Fermi level 11,12 . RKKY interaction on 
the graphene and graphcnclikc systems at half-filling has 
been constantly studied 11,12 , and the RKKY coupling on 
zigzag edge is reported to decay exponentially due to the 
localized zero-energy states 13 , while in graphene the cou- 
pling is proportional to 1/i? 3 . 14 

In this paper, we use the quantum Monte Carlo (QMC) 
method based on the Hirsch-Fyc algorithm 15 to study the 
indirect coupling between magnetic adatoms placed on 
the zigzag graphene edge. We calculate various thermo- 
dynamic quantities with respect to the chemical poten- 
tial, since graphene is a system whose carrier density can 
be tuned by the electric field effect 16 . For T k < T < D, 
where Tk is the Kondo temperature and D is the band 



width, the spin-spin correlation between the two mag- 
netic impurities are expected to be 17 
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where a Z i is the spin on the i-th impurity and /3 is the 
inverse temperature. Therefore our study on the spin- 
spin correlation between magnetic impurities can provide 
insights into the behavior of the RKKY interactions while 



II. RESULTS 

In our study, we choose zigzag edged ribbon with 10 
zigzag chains, where the density of states on the edge is 
comparable to the van Hove singularity at E = It in bulk 
graphene 5 . 

We fix the location of one impurity and change the 
location of the other impurity along the zigzag chain to 
study the behavior of the indirect interactions. Shown in 
Fig. 1 are the values of magnet moment squared, which 
is defined as 



m d = (i n dt - ridi) 2 ), 



(2) 



with nda being the carrier density on the impurity site 
and a being the spin index. When the two impurities are 
both located on the edge sites belonging to sublattice 
A, as shown in Fig. 1 (a), the local moment is strongly 
suppressed near the zero energy. Due to the indirect ex- 
change between impurities, when two impurities are lo- 
cated on Aq and A\ sites, as shown in Fig. 1 (c), the 
values of local moment is less suppressed. As the dis- 
placement between the two impurities is increased, the 
effect of indirect coupling decreases such that the local 
moment is more strongly quenched. The results for two 
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FIG. 3. (Color online). The weight of the different compo- 
nents of the wave function of d electrons on Aq and Ai. Here 
we use U = OM, V = 0.5t, /3 = 32i _1 and e d = -U/2. 




FIG. 1. (Color online). The local moment squared m\ as a 
function of /Lt. (a) Two impurities both located on the edge 
sites belong to sublattice A, (b) two impurities on different 
sublattices. Ao — Ai denotes the configuration that two impu- 
rities are located on Ao and A\ sites. Here we use U = 0.8t, 
V = 0.5t, /3 — 32t _1 and 6d = —U/2, (c) the schematic struc- 
ture of zigzag chain. One adatom is fixed on the red circled 
site Ao and the other one is moved along the zigzag chain. 
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FIG. 2. (Color online). The spin-spin correlation between 
two magnetic impurities on the zigzag edge as a function of 
li. One impurity fixed on the site Ao, the other impurity 
locatedon the site A\ (a), and Bo (b), respectively. Here we 
use = 32f _1 and e d = -U/2. 



impurities located on two different sublattices are shown 
in Fig. 1 (b). The effect of the indirect coupling is much 
weaker than the cases shown in Fig. 1 (a), and we can 
observe that the local moment on the A sites are not af- 
fected by the RKKY interaction, while that on the B site 
is slightly suppressed near the zero energy. The results 
for site configurations with larger distance are not shown 
because the effect of the indirect coupling decays so fast 
that the results would be mostly the same. 

Shown in Fig. 2 are the spin-spin correlation between 
two magnetic impurities with one impurity fixed on the 
site Aq and the other impurity located on the site A\, 
Bq, respectively. Various values of parameters as U = 
0.8*, 1.2* and V = 0.5*, 0.75* are tested for (3 = 32* _1 
and e<j = —U/2. In Fig. 2 (a), we can find that when 
fi = 0, the spin-spin correlation is ferromagnetic, which 
is consistent with the results of the RKKY interaction 



in the bipartite lattice 11,12 for the undoped case. As 
/i is lowered from zero-energy, the FM spin-spin corre- 
lation is first decreased and then at approximately at 
fi = — 0.15*, it vanishes for all the parameters we exam- 
ined. If we continue to decrease the chemical potential, 
the FM correlation turns to AMF and the magnitude 
grows as /i is decreased further. Hence, we can observe 
a transition from FM to AFM correlation, and the AFM 
correlation becomes stronger when we switch down the 
chemical potential. The AFM correlation is the strongest 
at (i ~ —0.3*, and the decrease in the value of /i will 
lead to decrease of the strength of AFM correlation. Fi- 
nally at /i ~ — 0.6t, AFM correlation vanishes. Near the 
zero-energy, the spin-spin correlations for U = 0.8* and 
V = 0.75* is more strongly quenched compared with the 
case U = 0.8* and V = 0.5*. As the hybridization V 
increases, the screening effect of the conduction electron 
would be enhanced, which will lead to much stronger sup- 
pression on the local moment formation. On the other 
hand, for the same values of hybridization V = 0.5*, the 
spin-spin correlation for U = 1.2* is larger than that of 
U = 0.8*, because larger repulsive U would give rise to 
single occupation on the impurity site. 

If the second impurity is located on the site Bq, as is 
shown in Fig. 2(b), the spin-spin correlation at undoped 
ribbon is antifcrromagnctic. However, as we switch fi be- 
low the zero-energy, the AFM correlation is enhanced, 
and it has a maximum magnitude at fj, ~ —0.25*. If we 
further decrease the chemical potential, the AFM corre- 
lation is decreased either, and finally at fi ~ —0.4*, it 
vanishes and the RKKY coupling changes to a FM one. 

The QMC method naturally returns the imaginary 
time Green's functions of the impurity sites, by which we 
can simply obtain the correlation functions such as (rii a ), 
(niaTijo), (n ia n' ilT n 3a ) and (n^n'^n^n'^). The two im- 
purity sites have 16 different spin configurations, and the 
sum of the weight of all these configurations is defined to 
be unity. All of the configurations would contribute to 
the spin-spin correlation, but our simulation reveals that 
the configurations | t, \) and | t, t) arc dominant near 
the zero-energy. 

In Fig. 3, we show the weight of different components 
of the d-electron wave function. As is shown in Fig. 3 
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FIG. 4. (Color online). The spin-spin correlation of the two 
magnetic impurities located on A-A and A-B sites in bulk 
graphene. We fix the location of one impurity on a site belongs 
to the sublattice A. A-A case means the other impurity is on 
the next nearest neighbor(belongs to sublattice A), and A-B 
denotes the other impurity is located on the nearest neigh- 
bor(belongs to sublattice B). Here we use U = 0.8t, V = 0.5i, 
(3 = 32i . (a) Two adatoms are located on the edge of 
ZGNR, (b) two adatoms are located on the bulk graphene. 



(a), the weight of | f, t) type dominates when \i = 0. As 
the chemical potential is lowered, we find that the weight 
of | tj t) type decreases and that of | f, l) is largely en- 
hanced. If we further lower the chemical potential, the 
contribution from the | tj 0) and the vacuum states in- 
creases due to the decrease in charge density. In Fig. 
3 (b), we notice the similar behavior. As \i is lowered, 
the weight of the | f, f) type wave function is drastically 
enhanced, while those of the | f, \.) is decreased, and the 
two curves cross over each other at fj, f» — OAt. Hence, the 
studies on the weight of different components of the wave 
function are in agreement with those of spin-spin correla- 
tions. As we switch down the chemical potential, parallel 
spin configuration becomes more and more favored by the 
two impurities. As \x is further lowered, the charges move 
out from the impurity such that the wave functions for 
the single occupancy and vacuum states dominate and 
hence, the AFM correlation is gradually decreased. 

Shown in Fig. 4 are the comparison of spin-spin cor- 
relation between two magnetic impurities located in the 
zigzag ribbon and in the bulk graphene system. In un- 
doped pure graphene as shown in Fig. 4 (b), when the 
two adatoms are located on the same sublattice, we can 
still see the FM correlation, but the strength is much 
smaller than that in Fig. 4 (a). As we lower the values 
of fi, the FM correlation is decreased, and at \i ~ — 0.2t, 
the correlation vanishes. For the A-B case, we find the 
AFM correlation is the strongest at fi = 0, and the val- 
ues of the AFM correlation is decreased as we switch the 
chemical potential away from the zero-energy. The be- 
havior of the spin-spin correlation in bulk graphene is 
completely different from that in zigzag edge of graphene 



ribbon. As we switch down fi, there exists an interval of /i 
that the correlation increases for decreasing fx. In zigzag 
graphene ribbon, due to the zero-energy states, the local 
moment formation at half-filling is strongly quenched, 
while in bulk graphene, the vanishing DOS at the Dirac 
cone gives rise to a well developed local moment. 

III. SUMMARY 

We examine the indirect coupling between two mag- 
netic impurities located on the zigzag edge of graphene 
ribbon, and compare the results with those obtained in 
bulk graphene sheet. Due to the zero-energy localized 
states on the zigzag edge, the local moment formation 
is strongly quenched. Interestingly, when the two im- 
purities are located on the zigzag edge, we find an en- 
hancement in AFM correlation with decreasing fi, which 
dose not exist in bulk graphene. The spin-spin corre- 
lation between two edge sites (A-A) is much enhanced 
when comparing with its bulk graphene counterpart, in- 
dicating the increase in indirect coupling caused by the 
zero-energy localized states. 
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